Fractals and Scars on a Compact Octagon 



Janna Levin and John D. Barrow 
DAMTP, Cambridge University, Silver St., Cambridge CBS 9EW 

and 

Astronomy Centre, University of Sussex, Brighton BNl 9QJ, U.K. 



A finite universe naturally supports chaotic classical motion. An ordered fractal emerges from the 
chaotic dynamics which we characterize in full for a compact 2-dimensional octagon. In the classical 
to quantum transition, the underlying fractal can persist in the form of scars, ridges of enhanced 
amplitude in the semiclassical wave function. Although the scarring is weak on the octagon, we 
suggest possible subtle implications of fractals and scars in a finite universe. 
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Chaos is often avoided in fundamental physics and cos- 
mology. A recent surge of interest in a finite universe has 
highlighted a well known chaotic system, a compact neg- 
atively curved space ||^. Most efforts in cosmology have 
struggled to circumvent the inherent chaos. Instead of 
avoiding chaos, we investigate the structure of the chaotic 
flows on a compact space. Since we are able to deduce 
analytically many of the chaotic properties, we concen- 
trate on the compact 2d octagon and study the classical 
motion and the emergent fractal properties. Poincare re- 
alized that the periodic orbits, though seemingly special, 
define the entire structure of a dynamical system. With 
the advent of chaos, periodic orbits grow exponentially 
in number as a function of their length. The invariant set 
of periodic trajectories packs itself into a fractal to ac- 
commodate the dense proliferation of orbits. We derive 
a very precise counting of the periodic orbits, determine 
the topological entropy, derive a recurrence relation for 
all of the fixed points, and analytically compute a spec- 
trum of fractal dimensions. We then numerically confirm 
our analytic results by generating the fixed points and 
measuring the dimensions of the set. 

The fractals may not be effaced in the classical to quan- 
tum transition. A phenomenon known as scarring of the 
quantum wave function has been witnessed in systems 
which are classically chaotic |^ . The scars refer to ridges 
along which the amplitude of the semiclassical wave func- 
tion is enhanced. We argue that the scars can be inter- 
preted as the smooth remnant of the underlying classical 
fractal. We relate the closed orbits on the 2d space to 
the occurrence of scars. The compact octagon is so thor- 
oughly chaotically mixed that the scarring is minimal, 
consistent with the fractal properties as discussed. In a 
finite model of the universe, the faint scars could seed a 
weblike distribution of matter threaded around the space. 

Geodesic motion on a compact 2d pseudosphere is com- 
pletely chaotic. On the Poincare unit disc the metric can 
be written 



-df]'^ + 4(1 - r2)-2 (^^2 ^ 
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Geodesies are semi-circles which are orthogonal to the 
boundary at r = 1. We consider null geodesies for which 
the coordinates (r, 0) have the conjugate momenta 



There are two conserved quantities. Since IT,^ 
is a conserved angular momentum 
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and the kinetic energy is conserved 
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Gondition (Q) is the same as the restriction that photons 
travel at unit speed. 




FIG. 1. The symmetric octagon drawn in the Poincare unit 
disc. The octagon forms a two-handled surface when the op- 
posite sides are glued together. Each edge corresponds to a 
generator Qj with j = —4, .., 3 and are numbered accordingly. 



We consider the simplest compact surface of negative 
curvature with genus 2. The fundamental domain is a 
symmetric regular octagon with the opposite edges iden- 
tified in pairs as in fig. |l|. The identifications are effected 
by a discrete subgroup of the Lorentz group. The hori- 
zontal identifying boost can be represented in coordinates 
X — r cos (f> and y — r sin (j) as 

_ / coshTi/2 sinhTi/2\ 
^° " i^sinhTi/2 coshri/2 j 

with cosh Ti/2 = 1 + \/2, while the other four are related 
by a rotation of go through the angle j7r/4. The four 
generators are related to their respective inverses by a 
rotation through tt. We denote all the generators and 
their inverses compactly as gj with j = —4, ..,3 so that 
the go5-4 = 53.9-1 = 929-2 = 9i9-3 = I- These genera- 
tors obey the one relation 



9o9-3929-i9-m9-293 
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A geodesic is uniquely specified by two parameters 
{6, L) where is the angular coordinate of the geodesic as 
it intersects the boundary in the past. Balasz and Voros 
Q developed a reentry map which takes initial values for 
(6*, L) through the appropriate edge and back into the 
fundamental domain. The map which iterates all values 
through edge of fig. |] is 



tan(6'„+i/2) = e^^^ tan(6l„/2) 
f 1 = (coshri + cos^n sinhri)L„ 



(7) 



sinh Ti sin 6„ 



The reentry map for orbits exiting any of the 8 edges is 
effected by shifting 0^9^ juj^ for j = -4, ...3. The 
map is Hamiltonian and area preserving. 

A trajectory is a point in the plane (0, L) and those 
trajectories lying within the intersection of the two dark- 
ened separating curves of fig. || exit edge 0. A separating 
curve is 



sin(6l - 7r/8) 



cos(6' — 7r/8) — cothr2 



(8) 
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FIG. 2. The action of the map on the exit-entry domains 
for geodesies. An orbit is completely specified by the angular 
momentum L and the angular coordinate ^(t = Q) = 9 a,t 
which the geodesies intersects the boundary in the past. All 
the orbits within the darkened exit domain above exit edge 
of fig. 0. The map takes this exit domain into the darkened 
vertical entry domain. 

with cothT2 = 2-5/'*(l + \/2). The other 7 exit domains 
are shifted by 7r/4 as represented in fig. ^j. All eight 
curves, Lj, are obtained from eqn. (||) by shifting 9 
9 + j7r/4 with j = —4, ..0, ..3. The map acts very much 
hke an eightfold Baker's map. The 8 sections exit their 
respective edges and reenter as the equal-area, nearly- 
vertical strips. Each of the 8 exit strips overlaps with 
only 7 of the reentry strips. In other words, the points can 
reenter any of the other edges except the one just exited. 
The original area is chopped into 8x7 smaller areas, each 
of which continues to flow through the space and gets 
resegmented into 8x7x7 new areas. The total number 
of pieces at order n is then 8 x T^~^ . Covering the area 
with boxes of size e = (1/7)"/^ gives the box-counting 
dimension Dq = lim„^oo ln(7""i)/ ln(7"/2) = 2. The 
fractal flUs the allowed area. 

The number of periodic orbits is slightly smaller than 
the number of areas. The gj construct a symbolic alpha- 
bet out of which to build n-letter words. Each repeated 
word corresponds to a periodic orbit. The number of 
word^ of order n, N{n), is the sum of the number of 
words which begin and end on the same letter, 0{n), 
and the number of words which begin and end on differ- 
ent letters, A(n), 



^The number of words actually counts fixed points and not 
periodic orbits. For instance, the words 171172 and g2gi describe 
two different fixed points along the same periodic orbit. These 
both contribute to N{n). The number of periodic orbits is 
slightly less than this, ~ N{n)/n. 



iV(n) = 0(n) -t- A(n). (9) 

A word cannot begin and end on inverse letters: a re- 
peated n-letter word of the form gigj---gkgi~^ is equivalent 
to the repeated {n — l)-letter word gj...gk. An n-letter 
word can be built from every 0(n — 1) by attaching any 
of the letters in the alphabet except the one which is in- 
verse to the end letter. So, there are 70(n — 1) of these. 
More n-letter words can be built from every A(n — 1) 
by attaching any of the letters in the alphabet except for 
the two letters which are either the inverse of the end 
letter or the inverse of the beginning letter. So, there 
are 6A(n — 1) of these. Although (n — l)-letter words of 
the form gi....g^^ do not contribute to N{n~ 1), they do 
contribute to A'^(n) (for n > 3) by adding any letter other 
than gi or g~^. The number of words of order (n — 1) 
that begin and end on inverse letters is 

m— n — 2 

/(n~l)= J2 A(m). (10) 

m=2 

It follows that 

N{n) = 70{n - 1) + 6A(n - 1) + 6/(n - 1). (11) 

Eqn. ( prj ) can be simplified by noting that 0{n) ~ N{n — 
1). We can combine these relations to reexpress (|ll|) for 
n > 3 as 

N{n) = 77V(n - 2) -t- 67V(n - 1) - 6A^(1). (12) 

There is additional pruning at higher orders due to the 
one relation eqn. (^. However, this pruning is extremely 
sparse and affects the counting very little. The number 
in (|l^ ) is actually slightly smaller than N{ri) — 7" but 
tends toward this behaviour as n ^ 00. To see this, 
let iV ~ &" in the limit of large n and solve for b to 
obtain N ^ 7". The topological entropy is then, Ht = 
lim„_oo (1/"-) In = In 7. 

There is a spectrum of weighted dimensions which 
measures the distribution of frequencies and the distri- 
bution of the points in space. The spectrum is defined 
by covering the set with boxes of size e on a side and 
taking the limit as as e — > 0, 

i?, = ^limi^ilf^. (13) 

q — I e^o Ine 

with Pi the weight assigned to the ith box. If all closed 
orbits are equally weighted then pi — 1/N{e) and the 
fractal is self-afiine with the Dq equal. For a multifractal, 
some of the orbits are more popular than others and at 
least some of the Dq are unequal. 

Following the reasoning of ref. Q for the Baker's map, 
we can use the areas of fig. || to analytically determine the 
spectrum of dimensions for the map. A given exit section 
of area A overlaps with the reentry strip into 7 overlap 
sections of area Aj . The map then takes these 7 areas and 
lays them down in nearly vertical strips preserving the 
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area of each. The natural measure on the vertical strip 
is then the relative area fj = Aj/A. The vertical strips 
will have nearly the same height and differ only in relative 
width, so the relative width is also given by fj . There is a 
simple scaling symmetry: if any of the vertical strips were 
to be multiplied by the width of that strip, the features of 
the entire set are reproduced. This kind of reproduction 
on smaller and smaller scales is precisely what makes the 
set fractal. We can exploit the scaling to deduce the 
fractal dimension j7|. In this vertical approximation the 
dimension is fractal only in the direction and Dq ~ 

1 + D^. Write the weighted sum of @, / = Yf'"^ Pt^ 
as /(e) = X]J=i -0 with natural measure fj on the jth 
strip. As we make e smaller by t/ fj, we reproduce the 
structure of the entire set. The scaling can be expressed 



A TST TjT TiT TST rjr TT^ 
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I,(6) = ///(e//,). 
From the definition (13), /(e) ~ gC^-i)^,^ -^^g obtain 



(14) 



3 = 1 



1. 



(15) 



Since the natural measures are normalized so that 
J2j fj = 1 it must be that the exponent q+ (1 — q)Dq = 1 
and it follows that 



(16) 



for all q. The fractal is self-similar and fills an area. It 
should be noted that this dimension is computed by tak- 
ing the entire area filled with periodic trajectories of all 
lengths and kneeding those areas into an iterate of par- 
allel sections, analogous to a Hamiltonian Baker's map. 
We can find a slightly different dimension which is sen- 
sitive to the number of windings of the orbit around the 
space. Numerically we generate the fixed points order by 
order instead so the area is initially empty and fills with 
fixed points at each iteration. We must get the same box- 
counting dimension for this set, that is, the set fills the 
allowed area, but the spectrum of weighted dimensions 
is different as shown below. 

To obtain the dimensions numerically we isolate the 
periodic orbits order by order analytically. The peri- 
odic orbits correspond to fixed points of the map. The 
points are then numerically covered with boxes in the 
{0, L) plane to explicitly determine the dimension. Let 
y — exp(i6'„+i), x = exp(«^„), and bj — e-np {inj /A). 
With a = tanh(ri/2) = ^2(1 + \/2)-i/2 the 6 map eqn. 
(^ can be written as 



y = 



abjX 



(17) 



for all points (0, L) in the jth diagonal strip exiting the 
jth edge. All n = I fixed points satisfy b']y'^ — 1 which 



FIG. 3. Each point corresponds to a periodic orbit. All 
fixed points up to order n = 5 are shown. There are 19, 624 
in all. 



implies 6', — ^ and i = 0. These correspond to the 4 
short period orbits and their inverses shown in fig. [l]. 

The fixed points at any order will satisfy an equation of 
the form y = {Ay + B)/{Cy -\- D) and the periodic orbits 
are solutions to a quadratic equation. The fixed points 
are found by using the following recurrence relations: 

^(n) ^ +05-1 n("-i) 

^ £,("-1) + abj(^)B^"-^'> . (18) 



At order n = 5, there are 



n— 5 



= 19, 624 periodic 
These are shown 



points exactly as predicted by eqn. 
in fig. |. 

Numerically we count the number of boxes needed to 
cover the set as e ^ 0. We measure the box-counting 
dimension as 2, as predicted. To find the weighted di- 
mensions we let Pi equal the number of fixed points in 
a box and it is important to note that we count each 
time a given fixed point occurs, so short orbits are more 
heavily weighted. We find the information dimension is 
Di — 1.4 ± 0.1. The weight assigned to each box is sen- 
sitive to the number of windings around the space and 
therefore measures a different topological property than 
eqn. dl^). 

In the transition to quantum mechanics, remnants of 
chaos have been found in the form of smooth scars of 
enhanced amplitude along the classical periodic orbits 
We argue that the scars are the residue of the most 
frequently encountered orbits in the fractal set. 

The semiclassical wave function can be constructed as 
a superposition of modes ^'(x) = X]fc '?^fcV'fc(^') with ipk 
eigenvectors of the Laplace-Beltrami operator, V'^ipk = 
— k'^ijjk- The intuition that the quantum system blurs 
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any chaotic features derives from two important conjec- 
tures. The first conjecture asserts that the eigenvalues 
are given by random matrix theory for systems which 
are classicahy completely chaotic. As a consequence, the 
amplitudes are drawn from a Gaussian random ensem- 
ble with a flat spectrum and should be thoroughly mixed 
and featureless. The second conjecture, due to Berry 
implies that the eigenmodes will be concentrated in the 
regions of phase space which are visited by the classical 
trajectories given an infinite time. A completely chaotic 
region of phase space is covered by mixing trajectories 
and at first glance Berry's conjecture seems to imply a 
uniform distribution. 

While these arguments for uniformity are convincing, 
one might wonder what happened to the structure of the 
periodic orbits. It is no surprise that the periodic orbits, 
which play such a prominent role in the classical dynam- 
ics, should also play a defining role in the quantum sys- 
tem. Periodic loops have in fact been the backbone of a 
quantization prescription in quantum chaos Scars 
are the first evidence of quantum features related to these 
closed orbits. Heller discovered scars in the evolution of a 
wave-packet moving along a short orbit in a chaotic sta- 
dium [|| . He deduced that only certain eigenstates would 
contribute to the wave function as it traversed the orbit. 
Since the wave-packet was peaked about a single periodic 
orbit, the few contributing eigenstates must also show en- 
hanced intensity along that orbit. The fewer the number 
of eigenstates which contribute to the wave-packet, the 
larger the intensity of the scar along that loop. 

Rather than being inconsistent with Berry's conjec- 
ture, the scar can be interpreted as a consequence of it. 
Scars occur in the regions of space where classical trajec- 
tories spend the most time. A typical trajectory will near 
a periodic orbit and then, due to the instability, will veer 
away onto a different periodic orbit. In this way, typical 
paths can be built from segments of periodic orbits. The 
most frequently visited periodic orbits tend to scar the 
most substantially. For a fully chaotic system, the frac- 
tal set of closed orbits will fill a volume in phase space. 
This would lead to some uniformity in the distribution 
of classical orbits and hence their quantum counterparts, 
as Berry suggested. However, not all closed orbits are 
created equal. Some may be visited by typical aperiodic 
trajectories more frequently than others. 

The scarring phenomenon has been investigated pri- 
marily in terms of the dynamics of a wave-packet for test 
particle motion. Since the periodic orbits are unstable, 
a particle will depart from that trajectory by a factor 
of exp (—At) where A is the positive Lyapunov exponent. 
The frequency uj with which the wave-packet traces out 
the path must be sufhciently short, X/oj < 1, if the orbit 
is to be susceptible to scarring. For paths of length £, the 
frequency is w = 27r/£. A hyperbolic surface is specified 
by A = 1 in curvature units and the scarring condition is 
simply £ < 2tt [§J|]. Only the shortest orbits scar. The 
short orbits are also the ones visited most frequently. The 



paths shown in fig. ^ are all of length £ 3.06 and are 
all potential avenues for scars. 

On the compact octagon, the eigenmodes must obey 
precisely the same conditions as the periodic orbits: 
i-'kix) — ipkigjx). The eigenmodes must be built out 
of primitives which correspond to the primitive periodic 
geodesies. For compact hyperbolic spaces, a relation be- 
tween eigenmodes and closed loops is evident. Aurich and 
Steiner investigated scars in an asymmetric octagon in 
Refs 1^ . They confirmed that the weighted spectral den- 
sity of states shows bands related to the periodic recur- 
rence of the wave packet along the orbit. The eigenstates 
under the bands had de Broglie wavelength, A_b ~ 27r/fc, 
a multiple of the length of the orbit, leading to the band 
quantization condition k ~ for integer band quantum 
number n ^ . The fewer the number of states contribut- 
ing to a band, the more likely they were to show scars. 
Consistent with this observation, the authors found that 
very low energy eigenstates did scar short orbits. The 
low energy eigenstates had few orbits to choose from and 
therefore tended to show scars. (They found scars can 
appear as underdensities as well as overdensities |^ . ) For 
higher energy eigenstates however, they found that scars 
were absent. High- A; states have many orbits to choose 
from and so are less likely to scar along a given one. In 
fact, since Dq = 2, we know that as A: — > oo, the area 
is covered by closed orbits. So, even if a given high-fc 
eigenstate aligns with some closed orbits, it will be en- 
tirely consistent with a random statistical distribution. 

We should emphasize that a full description of the 
eigenmodes requires an infinite sum over orbits, or at 
least a very large and cautiously truncated sum. Still, we 
can guess the location of scars in physical space without 
this detailed description by exploiting the fractal connec- 
tion. The scarred features can be described qualitatively 
for any classically chaotic system without having to quan- 
tize the system: simply isolate the most frequently visited 
orbits in the fractal set of periodic orbits. These will be 
the site for quantum scars. This offers a simple indicator 
of the quantum remnants of chaos. 

A more strongly multifractal set would show more ex- 
treme scars than the weak scars on the octagon. But the 
compact space is interesting as a model of the universe 
and, like our universe, is suggestively marginal between 
homogeneity and inhomogeneity. If the octagon were a 
2d universe, created in a Big Bang, then fluctuations in 
the geometry could be expanded in terms of the eigen- 
modes with amplitudes drawn from a random Gaussian 
ensemble with a flat spectrum. The superposition of fluc- 
tuations would show weak scars along short orbits from 
the low-k modes but would appear increasingly random 
on smaller and smaller scales. 

Since galaxies have their origin in quantum fluctua- 
tions, the distribution of large-scale structure would re- 
flect the homogeneity on nearly all scales except the 
largest where the pattern of galaxies could carry memory 
of the quantum scars. On scales comparable to the area 
of the universe, the distribution of matter could appear 
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filamentary and wcb-likc. 

Our own 'id universe does sliow vast voids and huge 
filaments in the distribution of large-scale structure. 
Galaxies would trace the scars and chart the shortest 
loops around a finite universe. If our universe is finite 
and negatively curved, the cobweb of luminous matter 
might be a residue of primordial quantum scars. 
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